Dissipative dynamics of magnetic solitons in metals 
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Soliton dynamics in spin-textured metals generate electrical currents, which produce backaction 
through spin torques. We modify the Landau-Lifshitz-Gilbert equation and the corresponding soli- 
tonic equations of motion to include such higher-order texture effects. We also find a quasistatic 
equation for the induced electrochemical potential, which needs to be solved for self-consistently, in 
the incompressible limit. As an example, we consider the gyration of a vortex in a point-contact 
spin valve, and discuss modifications of orbit radius, frequency, and dissipation power. 
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The current-driven dynamics of magnetic solitons such 
as domain walls, spirals, and vortices has recently at- 
tracted much attention!^ It i s actively researched for po- 
tential applications in magnetic random-access memoryP 
data storage devices,^ and high-quality tunable mi- 
crowave emitters.^' Soliton dynamics might also be rele- 
vant for current-induced switching in single-domain mag- 
nets, where magnetic textures are often nucleated in in- 
termediate stages during the switching process. Spin- 
valve structures provide an experimentally accessible 
arena for studying the current or magnetic-field driven 
dynamics of solitons or multi-solitonic statesPlt is there- 
fore desirable to have a precise theoretical understanding 
of the magnetization dynamics in spin-textured metallic 
magnets. 

In the literature, current-driven magnetization dynam- 
ics are usually described by the Landau-Lifshitz-Gilbert 
(LLG) equation with spin torques, or the corresponding 
equations of motion for collective coordinates describ- 
ing solitons. The electrical currents are conventionally 
treated as an external source and magnetic spin tex- 
ture is only taken into account through the exchange 
energy and the associated effective field. However, it 
is now established that magnetic texture dynamics gen- 
erate charge currents,^ which produce feedback on the 
magnetization via spin torquesP The three-dimensional 
coupled spin-charge hydrodynamic equations, including 
dissipation and thermal noise, were derived in Ref. [5]. 

The purpose of this paper is to theoretically exam- 
ine the role of such nonequilibrium spin-torque effects on 
the dynamics of the nanoscale spin-texture solitons. We 
find, in particular, that the magnitude of the corrections 
to the magnetic damping can be significant compared 
to the usual Gilbert damping in conventional transition- 
metal based materials. The additional terms in the LLG 
equation are nonlocal and anisotropic and cannot in gen- 
eral be absorbed into existing parameters. In this regard, 
they fundamentally modify the dynamics. Furthermore, 
because they are sensitive to the details of the magnetic 
order, they may be used to probe the local geometry and 
global topology of the magnetic texture. 

The LLG equation governing long-wavelength spin- 
texture dynamics in metals is given bjM 

S(rh + am x m) +m x H = V(l +/3mx)(j ■ V)m, (1) 



where we denote the local orientation of the spin-density 
field by the unit-vector field m(r,t), S is the magnitude 
of the saturated spin density, and j is the charge-current 
density. The magnetization vector is given by M = -fSm. 
In terms of the electron charge — e, V = ph/2e is the 
charge-to-spin conversion factor, where p is a dimen- 
sionless material-dependent spin-polarization parameter, 
which is typically a number of order 1. Here, H stands 
for the effective field defined by the functional derivative 
H = 5 m V, where V[m] is the free energy of the ferro- 
magnet. In a planar geometry, the in-plane current j in 
Eq. (JlJ has to be solved for self-consistently, taking into 
account spin-texture generated electromotive forces!^* 
Any additional spin torques injected vertically by exter- 
nally polarized spin currents will be included as contri- 
butions H st to the effective field. Equation (JlJ has been 
used to analyze the current-driven dynamics of magnetic 
vortices, spirals, and domain wallsP^ 

To illustrate the essential physics in the simplest set- 
ting, we start by considering the highly-compressible 
limit, such that the spin-charge dynamics may be de- 
coupled, resulting in a modified LLG equation (with the 
opposite limit discussed later) P^ESl 



{l + mx K) + mx (a + T) m + mxH/S = 0. (2) 



The current-texture interaction generated a correction to 
Gilbert damping given by 

F = — [(m x Vjin) <g> (m x Vim) - B 2 V l m ® V t ml , 
pS L J 

as well as a renormalization of the gyromagnetic tensor: 
++ BV 2 

K = — — [(m x Vim) ® V,m — V^m ® (m x V^m) . 
pS 

Here, p is the resistivity, the summation of the repeated 
spatial indices (i = x, y, z) is implied, and both tensors 
are given to quadratic order in texture gradients.. The 
total damping tensor determines the dissipation power: 

P= - [d 3 rH rh = S [d 3 rth- (a + T) m. (3) 



Let us estimate the size of the correction to the usual 
Gilbert damping using typical values for transition met- 
als: resistivity p ~ 100 Q nm, damping parameters 
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a,j3 ~ 10 2 , gyromagnetic ratio 7 « — 1.8 x 10 11 rad/sT, 
magnetization M ~ 10 6 A/m, exchange constant A xc ~ 
10~ n J/m, exchange length A = ^/ A xc / /iqAI 2 ~ 3 nm, 
and polarization factor p ~ 1. We find that numeri- 
cally, |T| - ^V 2 jpM\ 2 ~ a. Therefore, for textures 
with length scales on the order of A, the nonlocal damp- 
ing correction should be treated on equal footing with 
the usual Gilbert damping, the precise magnitude of the 
former depending on the material and texture. 

Numerical considerations aside, the nonlocal tex- 
ture effects give the LLG equation qualitatively dif- 
ferent structure. Solving Eq. ([2| directly is generally 
a formidable task that requires numerical integration. 
However, one is often interested in the dynamics of mag- 
netic solitons which are particle-like objects whose mo- 
tion may be captured by collective coordinates. One can 
then approximate Eq. ([!]) b y an equation of motion for 
these collective coordinates! 13 * 14 ! In the following, we add 
the nonlocal texture effects to the latter. We thus con- 
sider magnetic textures m(q(f),r) whose dynamics are 
parametrized by generalized coordinates q(t) . For exam- 
ple, they may be the location of topological defects such 
as vortices and domain walls in rectangular or curvilinear 
coordinates, or the out-of-plane angles and widths of the 
Neel domain walls, or the direction and pitch of spirals. 
The resulting equations of motion for the generalized co- 
ordinates will be generally nonlinear. 

Denoting the spherical angles of the spin density field 
m(q, r) by (6,ip), a general magnetic texture is charac- 
terized by the following second-rank tensor densities: 

bij = m • (d qi m x d qj m) = sm9(d qi 9d qj (f - d qj 9d qi ip) , 
dij = d qz m- d q . m = d q% 9d qj 9 + sin 2 9 d qi <pd qj ip , (4) 

Since b is the Jacobian d(— cos9, ip)/d(qi, qj) of the map- 
ping (qi,qj) —> (— cos#, (p), if the mapping is nonsingu- 
lar, its integral Gij = iwQij can be interpreted in terms 
of the number of times Qij the magnetization maps the 
{.Qiilj) plane onto the sphere. If the boundary condi- 
tions dictate to take discrete (e.g., integer) values, 
it becomes a topological invariant, which characterizes 
topological sectors of the spin textured If we take the 
partial derivatives in Eqs. Q with respect to the ordi- 
nary spatial coordinates r, b^/An is called the topolog- 
ical or skyrmion charge density and Qij the topological 
charge. Geometrically, b^ are the oriented surface ele- 
ments spanned by the gradients of the mapping from the 
collective coordinates to the magnetization unit sphere, 
q — > m. The symmetric tensor d determines the rate at 
which the local magnetization changes as a function of 
the generalized velocities, by m 2 = dijqiqj, which is pro- 
portional to the usual Gilbert dissipation power density. 

Consider first the case without damping, i.e., a — 
(3 = in Eq. ([!]). This Landau-Lifshitz equation with 
reactive spin torque follows from the variational prin- 
ciple, 1x1 x S m S — 0, with a magnetization action 
S[m] — J dtd 3 r£. 1 If instead of varying over the space 
of all possible magnetization fields, we restrict our vari- 



ation to the subspace parametrized by the generalized 
coordinates q, constraining the dynamics accordingly, 
m = qid qi m, we find an equation of motion for q given 
by J d 3 rd q .m ■ (m x [Eq. Q]). We use the same equa- 
tion to derive the equation of motion for the generalized 
coordinates from Eq. (J2|, including the dissipative terms: 

[G(q) + G'(q)]q + [£>(q) + £>'(q)]q = F(q) , (5) 

where the well-known LLG tensors are given by 

F i = -\j d 3 rd qi m ■ H st - ±d qt J d 3 rV , 

G = J d 3 r"b, D = aJ d 3 rd. (6) 

The spin torque corrections to the gyrotropic and damp- 
ing tensor are given by G^-(q) = J d 3 rd qi m ■ K ■ d qj m, 

D^-(q) = J d 3 rd qi m-T -dq^m. The additional dissipation 
power P' from the correction to the damping tensor may 
be expressed as P' — D'^qiqj. The validity of Eq. ^ 
depends on the accuracy with which one parametrizes 
the dynamical magnetic texture with collective coordi- 
nates. Typically, there will be a hierarchy of hard and 
soft modes, determined by their respective relaxation 
timesP^l The softest mode is usually the rigid transla- 
tion of the entire texture. For simplicity, this is the only 
degree of freedom we will retain, letting R denote the 
position of the defect center. 

Within each topological sector, solitonic configurations 
minimize the free energy U — J d 3 rV/ S, which includes 
both exchange and nonlocal magnetostatic interactions 
sensitive to the shape of the ferromagnet. For simplicity, 
we assume that the size of solitons are much smaller than 
the system size, neglecting boundary effects, so that the 
internal energy does not depend on the position of the 
solitons. The only role that the internal energy will play 
is to determine the profile m(r — R(f)), which is assumed 
undistorted under spatial translations. Therefore, in the 
evaluation of the tensors in Eq. Q, we may substitute 
<9r = — d r . In the following, when evaluating the ten- 
sors Q, we will thus take the derivatives with respect 
to spatial coordinates: d qi = = —d Ti = — V,. The 
spin-torque corrections to Eq. (|6| thus become: 

{G',D'} = -jg J d 3 r{(3(bd + db), (b 2 + f3 2 d 2 )} . (7) 

To express Eq. ([5]) in a vector form, we define the gy- 
rovector G by Gi = e y Gjfc/2 (along with the analogous 
definition for b in terms of b) : 

R x (G + G') + (D + D')R = -d R U + F st (R) , (8) 

where Ff = J d 3 rV,m • H st /5. Equation ^ can be 
thought as describing a massless particle moving with 
friction in a magnetic field given by G (which is some- 
times called the Magnus force), and an effective potential 
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U(R) that includes exchange, anisotropy, and magneto- 
static energies, driven by an external force F st , which 
may or may not be conservative. 

As an example, we consider two-dimensional (2D) vor- 
tex dynamics on a circular disk of radius L of negligible 
thickness. The tensors below are to be understood as per 
unit length in thickness. We will denote the position r 
with respect to the center of the disk with polar coordi- 
nates (r, <fi), with the vortex center being at (R, ip). The 
vortex spin texture is described by 

V?(r') = tt/2 + 0' , 9(r') = 2 arctan(r'/A) atr'<A, (9) 

and 8(r' > A) = tt/2, where r' = r — R, whose polar coor- 
dinates are r' = |r — R| and 0'(R) = arg(r — R). Such a 
vortex profile crudely captures a compromise between the 
exchange and magnetostatic energies.^ In the coordinate 
system centered on the vortex, the texture gradients are 
given by V# = rd r 8, 'Vip — <fi/r, so that d rr — (d r 9) 2 , 
dcfcf = (sin6>/r) 2 , d r< p = 0, and b r( p — —{d r cos6)/r. The 
gyrovector density is b = b r( j,z, and the damping tensor 
density is d = d rr r ® f + d^<p ® <fi. Our vortex profile 
has the property that d r = sin^/r = f(x)/X inside the 
core, where x = r/X and f(x) = 2/(1 + x 2 ), so that 
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(//A) 2 at r< A; 



b r <f> = d rr = , d^ = 1/ (Ax) 2 at r > A . 



(10) 



In 2D, the gyrovector has only one nonzero component 
perpendicular to the plane, G = Gz, where G = G r $- 
Up to continuous deformations, it depends only on the 
vortex profile in the core, where our vortex is a soliton 
with topological charge Q = 1/2, so that G = 4ttQ = 2ir, 
as one can verify explicitly. Computing the spin-torque 
correction for the vortex profile ([9]), we find 



fit _ ft'P 

pSX2 



dx 4:77 X 



l + x 2 



(11) 



The integral evaluates to about 30. For the previously 
listed parameters, we find this 0(j3) correction to be 
small, \G'\/G ~ 10~ 4 , which can be safely neglected. 
Consider the damping tensor D in Eq. Q. Because 
we assume a rigid vortex which has rotational symme- 
try about its core, this tensor is diagonal: D = D = 
Tr[£>]/2 — a J d 2 r (Viin) 2 /2, which is proportional to 
the total exchange energy. Recalling that the vortex 
is far away from any boundaries, we find D — a[2n + 
7r ln(L/A)]. The 2n comes from the region inside the core, 
where the vortex |9| satisfies the minimum energy condi- 
tion Z7 xc /j4 xc = D/a = G. The logarithmic factor comes 
from outside the core and is generic to any 2D vortex. 
For the spin-torque correction to damping, we neglect 
the 0((3 2 ) contribution outside the core, so that 



D' 



V 



P SX 2 



dx 2ttx 



1 



(12) 



Taking disk of radius L ~ 10 3 A, we find that the damping 
is enhanced by D'/D - V 2 /apSX 2 - 1. 
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FIG. 1: (Color online) Left panel: Energy density (13 1 of the 



vortex texture m (red dotted arrows) due to the Oersted field 
H (green solid arrows), in the xy plane centered at the point 
contact, with unit length a. The arrows depict a crude model 
of the vortex as four domains separated by 7r/2 walls emanat- 
ing from the vortex center. On the gray scale, higher energy 
density corresponds with lighter shade. One clearly sees a 
region of high energy that scales linearly with the separation 
distance between the vortex and the point contact, resulting 
in a linear attractive potential. The inset shows a side-view 
schematic of the magnetic bilayer in the z direction. Right 
panel: A plot of the charge current due to each term in the 
spin-force field ( 19 1 in a region of 3A x 3A centered at the vor- 
tex core. The contribution collinear with the Magnus force 
is shown with thick red arrows. The viscous 0(/3) terms are 
indicated by blue arrows and have been magnified by setting 
/? ~ 1. The long black arrows indicate the direction of each 
term in Eq. (15 1, labeled by their respective couplings. 



To explore the consequences of our corrections, con- 
sider large-amplitude vortex dynamics in a point contact 
spin valve, driven by a dc current applied perpendicu- 
lar to the plane. 17 Electrons pass from a pinned ferro- 
magnetic layer with magnetization along the direction 
p, become spin-polarized along the same direction and 
exert a spin torque on a free layer with a 2D vortex. 
For simplicity, we assume the free layer is circular and 
p is uniform. Consider the potential energy of the vor- 
tex, LT(R) = - Jd 2 r H • M/S, due to the Oersted field 
H = —h(r)Icf), created by the charge current — Jz, where 
h(r) is the radial profile. Assuming the current is applied 
to a circular region of radius a, h(r) = p,o/2irr at r > a. 
The potential energy (per S) is thus given by 



U(R)= 1 I / drd(j>rh(r)smO(r')<t>- 4> 



(13) 



where <j> = z x f'. When the vortex is outside of the 
current distribution, the potential is approximately linear 
in See Fig. [l] We will therefore consider the linear 
potential U(R) ~ const + AIR, where on dimensional 
grounds, A ~ —70/10 ~ 10~ 3 m 3 /sA, taking a ~ 100 nm. 

In addition to the Oersted field, the current exerts a 
spin torque on the free-layer magnetization which drives 
vortex precession. We take a simple model, assuming 
that the transverse spin angular momentum of the polar- 
ized current is absorbed with an angle-independent effi- 
ciency a by the free layer. This spin torque can then be 
included in Eq. (pi by add ing t o the effective field a term 
H s t = (aVI/nar)p x mPES Inside the point contact, 
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Vim w — sin^Vi^f' (when the vortex core is outside), 
so that only the out-of-plane component of the polariza- 
tion p± = p • z, which creates an in-plane effective field 
along i x m = — sin 9r ', contributes to the resulting force 



to 2 .!? 



^rflin^Vpw-^^, (14) 



where ij> = i x R and R = RR. Here, we have taken 
the limit R 3> a, approximating the integrand with its 
value at the origin. Defining B — aVpj_/S, Eq. Q in 
the region R > a becomes 



Gz x R - L>R - AIR, - Biz x R/i? = , 



(15) 



where we absorb D' in D hereafter. It is natural to switch 
to complex notation, Z = X + iY = Re % ^ , and solve 



Eq. ( 15 ) in polar coordinates: 



R + iRtp 



T {AD — BG/R) + i(AG + BD/R) 

iRW — —1 — s -rr: 



D 2 + G 2 



(16) 



A circular orbit, R — 0, exists if sign(£>G) > 0, which 
for our vortex (G > 0) requires sign(i?) = sign(pj_) > 0. 
The orbit has radius Rq = BG/AD and frequency ip = 
—AI/GRo, and is stable only for I > corresponding 
to an attractive potential U(R). A twofold texture en- 
hancement of D would halve the orbit radius and double 
the frequency and dissipation power, P — D(AI/G) 2 . 

We now return to relax the assumption of high com- 
pressibility, which underlies derivation of Eq. ([2]). The 
local Ohm's law is given bj0 



P3i 



ViPVtin + m x Vim) • m + Vip/e . 



(17) 



where p is the drude resistivity, p(n) = n/K is the chem- 
ical potential, n is the nonequilibrium particle density 
and K is the compressibility. In Eq. (|2j), we decoupled 
the current-magnetization dynamics by taking p — > 
and substituting Eq. (17) into ([lj. For a finite com- 
pressibility, however, there is a particle-diffusion current 
ju = 'Vp/ep that gives an additional spin torque te> 
in Eq. ([!]) ~* This spin torque couples the magnetization 
dynamics to the charge density, which is determined by 
the continuity equation, en = V • j. 

We consider a typical scenario in which Coulombic 
electron-electron repulsion renders electric flows essen- 
tially incompressible, and solve these coupled equations 



in the n — > limit. The continuity equation then deter- 
mines the electrochemical potential according to 



V 2 Ai = "PeV-F, 



(18) 



where we wrote Eq. (17) asj = —VF/p+j£>. The texture- 
dependent field F = e + f (force per unit of "spin charge" 
Ve) is due to the fictitious electric field e and dissipative 
field f with components = m ■ (m x V,m) and — 
— /3m- Vim. A rigid soliton in motion generates e = Rx b 



and f = /3d ■ R. Formally, Eq. ( 18 ) is Poisson's equation 
for a charge distribution with polarization F, and one 
can apply standard Green's function methods to solve for 
p, once boundary conditions are specified^ Our vortex 
generates the force field 



F(r') 



(/3R + Rxz) (f(x')/\y 
/3(R ■$')<£' /(Xx')\ 



r' < A 
r' > A 



(19) 



See the right panel of Fig. [T] where the charge current 
generated by F is plotted. The divergence of this current 
induces a counterbalancing diffusion current ju, accord- 
ing to Eq. (18), which ensures that the total current is 



divergenceless. Neglecting the 0(/3) terms, the polariza- 
tion charge n po i is given by 



i P oi( r ') 



Ve 

47T 



V F 



Ve 
4ttA 3 



f'-(zxR)^(/ 2 ). (20) 



Since the charges vanish at the boundary, we obtain the 
Coulombic solution ^(r) = J d 3 r'n po i(r')/\r — r'|. 

In summary, we have extended the LLG phenomenol- 
ogy to include spin-texture effects stemming from dy- 
namically generated electric currents, and examined 
these effects in 2D vortex motion under an applied point- 
contact current. Our theory can be applied to multi- 
soliton dynamics, with or without applied currents and 
magnetic fields, in any dimensions. With a more detailed 
parametrization of the soliton profile, its equations of mo- 
tion may be rendered more realistic. The /3-like viscous 
coupling between current and magnetization dynamics 
can have giant enhancement)^! i n dilute magnetic semi- 
conductors due to the nonadiabatic spin-torque effects 
that are beyond the scope of Eq. (JlJ . This should moti- 
vate further study in the present context. 
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